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We present a lattice study of fluctuations and correlations among the conserved charges baryon 
number and strangeness in (2+l)-flavor QCD. The lattice calculations are based on a Taylor ex- 
pansion of the pressure. Results are presented at zero and nonzero density on lattices with four 
and six time slices, corresponding to a lattice spacing of a w 0.25 fm and a w 0.16 fm, respec- 
tively. The quark masses are almost physical, i.e. the light quark mass has been chosen to one 
tenth of the physical strange quark mass while the strange quark mass has been set to its physical 
value. We find that all analyzed fluctuations and correlations develop a peak at nonzero density, 
when they are treated in an expansion up to sixth order. Especially below the critical temperature 
(T t ) fluctuations and correlations increase drastically, whereas above T c they are rather unaffected. 
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1. Introduction 

Fluctuations of conserved charges, like baryon number, electric charge and strangeness are 
generally considered to be sensitive indicators for the structure of the thermal medium that is pro- 
duced in heavy ion collisions In fact, if at non- vanishing baryon number a critical point exists 
in the QCD phase diagram, this will be signaled by divergent fluctuations of e.g. the baryon number 
density [@]. 

We present here results from lattice calculations of baryon number and strangeness fluctuations 
in QCD with dynamical light and sttange quark degrees of freedom. The results are based on 
calculations with an improved staggered fermion action (p4-action) that strongly reduces lattice 
cut-off effects in bulk thermodynamics at high temperature. The values of the quark masses used 
in this calculation are almost physical; the strange quark mass, m s , is fixed to its physical value 
while the light up and down quark masses are taken to be degenerate and equal to m s /lO. This 
is about twice as large as the average up and down quark masses realized in nature. We obtained 
results from calculations performed with two different values of the lattice cut-off, corresponding 
to lattices with temporal extent N Z = A and 6. This allows us to judge the magnitude of systematic 
effects arising from discretization errors in our improved action calculations. The spatial volume 
has been chosen to be V^ 3 T =4, which insures that finite volume effects are small. 

At the QCD critical endpoint (CEP) the correlation length of the chiral critical mode a will 
diverge. Correspondingly, all kinds of fluctuations of conserved charges, which couple to the a- 
field, become large in the vicinity of, and diverge at the CEP. It has been argued in ^ |[] that 
quadratic variances of event-by-event observables such as particle abundances, particle ratios or 
mean transverse momenta will reflect these divergent fluctuations and are thus good experimental 
observables for the determination of the CEP. It will thus be very interesting to study baryon number 
and strangeness fluctuations as well as baryon number-strangeness correlations at nonzero density 
as they are related to the event-by-event fluctuations of the proton, the kaon and their ratio. The 
latter has been recently measured by the NA49 collaboration for central Pb-Pb collisions at five 
different SPS energies [§]. 



2. The Taylor expansion method 

Direct lattice calculations at nonzero baryon density are not possible by means of standard 
Monte Carlo methods. We follow here the Taylor expansion approach to finite density QCD on the 
lattice, as described in detail in [§,§]. Starting from an expansion of the logarithm of the QCD 
partition function, i.e. the pressure, one obtains 

^ = ^mz(y,r,^^) = Ic«f ) ; ' , (2.i) 



where 

B<s . . 1 d< dJ ]nZ(V,T,n B ,n s ) 



c 



tf(T) 



(2.2) 

Ms=Ks=0 



dn> B dn J s vr 3 

Here Hb and Us are the baryon and strangeness chemical potentials, respectively, which can be 
obtained as appropriate linear combinations of the quark chemical potentials. Note that we treat 
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Figure 1: Diagonal (upper row) and off-diagonal (lower row) expansion coefficients of the pressure of 2nd 
(left column), 4th (middle column) and 6th (right column) order. For diagonal coefficients up to the 4th 
order we compare results from two different lattice spacings (N% = 4 and N T = 6 lattices). Here we have 
normalized the data by its corresponding SB-values. 



the up and down quarks as degenerated, both in mass as well as net quark number density. Accord- 
ingly, isospin and electric charge chemical potential vanish. Due to charge conjugation symmetry 
only coefficients with i + j even are nonzero. In Fig. [I] we show the diagonal and off-diagonal 
coefficients of this expansion up to 6th order. For diagonal coefficients up to 4th order we show 
results that have been obtained for different lattice spacings, i.e. on N z = 4 and N x = 6 lattices. We 
note that results obtained on lattices with temporal extent N x = 6 are in good agreement with those 
obtained on the coarser N T = 4 lattice. The slight shift towards smaller temperatures, visible for the 
Nt = 6 data relative to the N x = 4 data, is consistent with findings for the equation of state, e.g. the 
trace anomaly (e — 3p)/T 4 , and also reflects the shift in the transition temperature observed when 
comparing the locations of the cusp in the chiral susceptibility [^]. 

The general pattern of the coefficients can be understood from the structure of the singular 
part of the free energy. In the chiral limit of our (2+l)-flavor simulations, i.e. with vanishing light 
quark masses, but finite and physical strange quark mass, we expect the QCD transition to be of 
second order [§J. The relevant scaling field which contains the baryon chemical potential will be 
the reduced temperature, for which we make the ansatz 



'=-tt +K yf> ■ <23) 

Accordingly two derivatives of the free energy with respect to fig (at Hb = 0) are similar to one 
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Figure 2: A fit to the expansion coefficients of the pressure below T c , inspired by the critical behavior of the 
free energy. At this lattice size and quark mass (N z = 4, m q = Q.\m s ) we have T c m 202 MeV. 



derivative with respect to the temperature. Therefore we obtain the following formulas, which can 
be used to fit the critical behavior of the Taylor expansion coefficients 

2!-45 ~ ^2A ± (2-a)K\t\ l - a + b 2 t + c 2 ; (2.4) 

4! -4,5 \2A±(2-a){l-a)K 2 \t\- a + b 4 t + c 4 ; (2.5) 

6! -eg ~ ±l20A±(2-a)(l-a){-a)K 3 \t\- l - a . (2.6) 

Here uppers signs are valid for t > 0, whereas lower signs are for t < 0. We see that the critical 
behavior is governed by the critical exponent a, which is small and negative for the universality 
classes of interest. For staggered fermions at finite lattice spacing, were the chiral symmetry is 
broken down to a £7(1) subgroup, we expect the relevant universality class to be that of the 3d-0(2) 
symmetric model, where we have a ~ —0.015. We thus expect Cg° s to develop a cusp, while Cg° s 

2 

will diverge in the chiral limit. c B ' s will be dominated by the regular part of the free energy. Hence, 
we took into account the leading order regular terms, indicated by the coefficients Z?2,^4,C2, and ca. 



in Eq. |2^ and Eq. |2.5| . For T < T c we show the fit results in Fig. g. Although the ansatz is strictly 
valid only in the chiral limit, the fit works quite well. It has already been observed previously that 
at this quark mass (m q = 0.1m. s ) the QCD thermodynamics enters the scaling regime of the critical 
behavior [Jsj] _ Currently the coefficient K, which determines the slope of the critical line as function 
of /ig is, however, not very well constrained. 



3. Baryon number and strangeness fluctuations 

We now construct baryon number and strangeness susceptibilities which measure fluctuations 
at nonzero chemical potentials. The corresponding expansions can be expressed in terms of the 



pressure coefficients from Eq. 2. 1 . We obtain for the baryon number fluctuations 

to = Mz^ = | Ca)(a _ I)<j!f (^, (3,) 
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Figure 3: Baryon number fluctuations (left), strangeness fluctuations (middle) and baryon number 
strangeness correlation coefficient (right) at nonzero baryon chemical potential. Data points are obtained 
by a truncated Taylor series up to the 4th order in Hb/T, while dashed lines indicate the 6th order results. 



and for the strangeness fluctuations 



Here (•) indicates the expectation value evaluated in the grand canonical ensemble in lattice sim- 
ulations. Nb and Ns are the net number of baryons and strange particles, respectively. In Fig. || 
(left and middle) we show results for the baryon number and strangeness fluctuations at zero and 



nonzero (J-b/T, obtained by truncating the series (3.1 and 3.2) after the 4th order. We have also 
indicated the 6th order results by dashed lines in order to give a feeling for the truncation errors. At 
Hb/T = 0, the generic form of the temperature dependence is a smooth crossover, while the leading 
order term in (hb/T) 2 has a peak at T c . It is thus clear that at sufficiently large values of {J-b/T the 
fluctuations will develop a peak. However, it is important to stay within the radius of convergence 
of the series, in order to keep the truncation errors small. We will discuss the radius of convergence 
in Sec. |]in more detail; here we find that for }Xb/T < 1.5 the truncation error is small. For this 
value of the chemical potential the baryon number fluctuations show already a small peak, whereas 
the strangeness fluctuations are only slightly enhanced. 

We note, that all fluctuations of modes that couple to the relevant chiral critical mode at the 
QCD critical point (the sigma-field) will diverge at the QCD critical endpoint (CEP). As light 
quarks carry also a baryon number, the baryon number fluctuations are expected to diverge at the 
CEP. 



4. The baryon number - strangeness correlation 

The baryon number-strangeness coefficient was introduced as a diagnostic tool for the effective 
degrees of freedom in the quark gluon plasma [^]. It is defined as 

r _ JNbN s )-(Nb)(N s ) 

Cbs=— 3 — - — -j . (4.1) 

(Nl)-(Ns) 2 
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Lattice calculations of this quantity in the (partially) quenched approximation [10] and 2 + 1 -flavor 
QCD [§] indicate that this quantity reaches unity for temperatures larger than l.5T c , whereas it is 
significantly smaller for T < T c . I.e., the quark flavors are uncorrected above T c as in the ideal 
Quark Gluon Plasma (QGP). This behavior indicates the transition from hadronic degrees of free- 
dom to that of a QGP. 



Similar to the expansion of baryon number fluctuations (Eq. 3.1) and strangeness fluctuations 



(Eq. p.2| ) we can also expand the baryon number- strangeness coefficient. Again we express its 
expansion coefficients in terms of the expansion coefficients of the pressure (Eq. |Q| ) and obtain 



with 



and 
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In Fig. || (right) we show the correlation Cbs at zero and nonzero density. We find that at 
nonzero baryon number density a peak is developing already at small values of Hb/T where the 
truncation errors are small (data points are the 4th order results, while the dashed lines indicate 
the 6th order results). At Hb/T = 1.0, the correlations at T c increase by almost 50%, and rise 
to a value which is about 20% larger than the Stefan-Boltzmann value. It will be interesting to 
see whether the increase in the baryon number strangeness coefficient will give rise to a peak in 
proton over kaon fluctuations obtained in Heavy Ion Collisions and whether these fluctuations can 
be used as an experimental signal for locating the critical point. Also below T c the correlations 
increase drastically with increasing whereas above T c the correlation coefficient Cbs remains 
rather unaffected by the increasing baryon density. 

5. The radius of convergence 



It is known that in the temperature region close to T c the convergence properties of the Taylor 



expansion of the pressure (Eq. 2. 1 ) and therefore also of the expansions of the baryon number and 



strangeness fluctuations (Eq. 3.1) and (Eq. 3.2) are poor. Here the expansion coefficients wildly 
fluctuate in sign and magnitude. Only for temperatures below the CEP, this might be better as in 
this case all expansion coefficients are positive. It has been argued that the radius of convergence 
can be used to determine the location of the CEP [|l2|]. One way to define the radius of convergence 



of the pressure series mpLs/T (Eq. |2T| ) is 



p = limp„ with Pn = H B n) /TW = Jc B n $/c B n 



B,S i B,S 

+2,0 • 



(5.1) 



The convergence radius can be estimated in a similar manner by the coefficients in the series of 
the baryon number fluctuations (Eq. |3~l| ). Each order p„ will, however, differ by a constant factor 
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Figure 4: Two estimates of the radius of convergence of the pressure series (full symbols) and the series of 
baryon number fluctuations (open symbols) as a function of temperature. p\ indicates our best estimate for 
the location of the critical endpoint. Dashed lines show the resonance gas limit of the different estimators. 



which goes to one in the limit n — > oo. We define a second estimator of the convergence radius, 
based on the ^-series as 

PnlXB) = J , n } n ^ l ] , Pn[p/T A ] ■ (5-2) 
y (n+\)(n + 2) 

The method to estimate the location of the CEP by the radius of convergence works in two steps: 

1. Find the largest temperature where all (available) expansion coefficients are positive. Only 
if all expansion coefficients are positive, the corresponding singularity which is limiting the 
convergence radius lies on the real axis and can be associated with the physical phase transi- 
tion. 



2. Estimate the radius of convergence at this temperature by using Eq. |5JJor Eq. p.2 . 



It is clear, that in practice we can not perform the limit n — ► oo but have to stop at some finite n. 
Eventually at this point different estimators of the radius of convergence will agree within errors. 

In Fig. |] we have summarized current results on the radius of convergence. The first non trivial 
approximation which does not explicitly include contributions from the gluon sector (contributes 
to c ' ), is given by P2, shown by the square data points. Here both contributing expansion coeffi- 
cients (<?2 'q and c 4 ' ) are positive for all values of the temperature, hence an approximation for the 
temperature (T^) of the CEP can not be deduced. At this order the estimates from the pressure se- 
ries (full symbols) and from the series of the baryon number fluctuations differ quite dramatically. 
The next higher approximation P4 of the CEP, which is currently our best approximation, is shown 
by the circle data point. We obtain for the temperature 0.96T c (/j,b = 0). Here the difference 

between the two estimators becomes less severe. For the convergence radius they span the range of 
p w (1.5 — 2.4). It is interesting to note that the two different estimators of the convergence radius 
seem to approach the n — > oo limit from different sites. However, the ^-estimator seems to be the 
more stable one. Also shown are the resonance gas limits of the two approximations, indicated by 
the dashed lines. 
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Figure 5: Baryon number fluctuations at T = 0.96T C w T L as a function of the baryon chemical potential 
(left). We compare different order of the truncated Taylor series vs two Pade approximants. The right panel 
shows the same but on a much larger scale. 



6. Taylor expansion vs. Pade approximation 

So far our analysis of observables at nonzero baryon density has been based on the Taylor 
expansion of the pressure (Eq. [O] ) around }Xb/T = 0. Now we want to discuss the Pade approxi- 
mation, which is usually considered to be a better approximation of the true function, even beyond 



the radius of convergence QllQ. One can construct different Pade approximants from the Taylor 
expansion coefficients. Pade approximants are rational functions which differ by the order of the 
polynomials in the numerator and denominator. For the baryon number fluctuations we obtain e.g. 

Ac 2fi c Afi +2 ( 12 (c^ -Sc^c 6 A (i^) 2 

^ t B,5 c B,5 ' z I lZl \ B,S J JC B.S C B,S \ T J 
Pade t 2 ' 2 ] = ^ 2 4,0 , 6,0,^2 (6-D 



In Fig. ||(left) we compare different orders of the truncated Taylor series (Eq. 3.1) with the Pade 



approximations [2,2] and [4,2], where the latter one has a fourth order polynomial in the numerator. 
We show the baryon number fluctuations at fixed, but varying baryon chemical potential. T 
0.96r c is currently our best approximation for the temperature of the critical endpoint (r CEP ) as 
it was estimated by the radius of convergence method. The most prominent feature of the Pade 
approximations is the appearance of a pole. The location of the poles are given by the roots of the 
denominator. From Eq. ^Tl] it is obvious that the pole in the [2,2] approximation exactly coincides 
with the P4 [xb] estimate of the convergence radius. 

Note that for the Taylor series of Xb which was truncated after the 6th order, as well as for 

8 

the [4,2] Pade approximant we need the coefficient c B S . In the following analysis we set the so far 
unknown coefficient to = Cg^/x. 

Given the coefficients Cg S - c 6 B ,0 s at T « 0.96T C , we vary the strength of c B '° s by the parameter 
x. We have fixed the parameter x by assuming that the asymptotic behavior of the [4,2] Pade 
approximant for large should be governed by the free gas result which is given by 
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Using this property, we obtain x 4.95 and the next approximation for the convergence radius 
(p6 = V*)> which we also shown in Fig. fright). Again we find that the Pole in the [4,2] Pade 
approximant for %b coincides exactly with the corresponding estimator for the convergence radius 

7. Summary and conclusions 

We have performed a comprehensive study of the Taylor expansion coefficients of the baryon 
number fluctuation, strangeness fluctuation as well as their correlations, based on the Taylor expan- 
sion coefficients of the pressure. We found that all these quantities will develop a peak at T « T c for 
nonzero baryon chemical potential, which is the dominant effect from the leading order coefficient 
in pLg. The sub-leading coefficient will in general lead to a shift of the peak, i.e. a /^-dependence 
of the peak position. Our current analysis takes into account the coefficients up to 6th order. On the 
highest order the errors are barely under control. As most observables at nonzero chemical poten- 
tial crucially depend on the relative strength of the coefficients, a more refined analysis which also 
takes into account the 8th order coefficients is highly desired. This will also help to confirm the 
structure of the Taylor expansion coefficients which can be understood in terms of an appropriate 
scaling ansatz of the free energy in the chiral limit. With a combined fit of all Taylor expansion 
coefficients it will in principle be possible to determine the /i B -dependent curvature of the critical 
temperature in the chiral limit. 

When comparing different orders of the truncated Taylor series we find small truncation errors 
below jUs/r<(l — 1.5), while for /J-b/T « (1.5 — 2.5) consecutive orders become compatible. The 
latter fact can be used to estimate the radius of convergence of the Taylor series, which will be 
connected with a physical singularity in the QCD phase diagram (critical endpoint) when all ex- 
pansion coefficients are positive. The temperature of the critical endpoint we currently approximate 
to T CEP « 0.96r c . 

We also constructed Pade approximants from the Taylor coefficients. Again we find good 
agreement between different approximants as well as the truncated Taylor series below /ig/r<(l — 

i.5)r c . 
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